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RECURSION FORMULAS FOR THE COEFFICIENTS OF THE 

f AND g SERIES 


If the xy-plane coincides with the plane of motion and if the positive x-axis 


is along the position vector r 0 at time 
in the form 

xy - yx 
xx + yy = r 0 x 

x ( t o) = x o = r o ’ 


t 0 , the two-body equations can be written 


y 0 r o > 

(i) 

o 

1 

<<* 

o 

o 

(2) 

y( t o) = y 0 - 0 - 



where /x is a constant and r is the magnitude of r. Zero subscripts designate 
values at time t 0 , and a dot over a symbol indicates a time derivative. 

Equation 1 expresses the constancy of angular momentum and Equation 2 is 
obtained by multiplying through by r the equation 


r - (/u.e sin0)/h 


- [(/ze sin0 o yh] cos (d -6 0 ) + [(/xe cos <9 0 )/h] sin [8 - d Q ) , 

where 6 is the true anomaly, e the eccentricity, and h the magnitude of the 
angular momentum. 

In the reference frame described above, the equation 

7 = f7 0 + gt 0 (3) 

reduces to 


x - fr 0 + gr 0 , 

y = gy 0 • 


(4) 

(5) 


1 


Equation 2 is not valid when y 0 = 0, i.e., for the rectilinear cases. Sub- 
stituting Equation 5, however, into Equations 1 and 2, 


xg - gx - r 


xx + bgg - r Q x + eg , 


b — v 2 — v 2 - r 2 

u y o v o r o 


b - /j./r 0 , 


Substitution of 


CP oo 


into Equations 6 and 7 with a 0 = r Q , aj = r 0 , b 0 = 0, and bj = 1 gives, for i > 1, 


l-l 

r o (i + 1) a i+i = a i a i + cb i " 2^ (j + X) ( a j+i a i-j +bb j+ i b i-j) - (!0) 


1- 1 

r o ( i + 1) b i+i = 2Z (i + !) ( a j+ i b i-j “ b j + i a i-j) • 


Taking advantage of symmetry in the summations, Equations 10 and 11 
become, for i > 2, 


r o a i+x 


k k 

+cb i)/( i + l)-2^ a jn -lW 2 +bb k«).< 12 ) 
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r o( i + 1 ) b i + i " 2[(i-J)-(i*l)](b j „ a ,,-ViV J ) . (13) 

3=0 

k = integral part of (A) (i - 2), 
q = fractional part of (A) (i - 2), 
a 2 = “M/2r 0 2 , b 2 = 0 

Having x(t) and g(t), f can be computed from Equation 4 and7(t) from 
Equation 3. Inversion of Equations 6 and 7 yields 

r 2 x = (f 0 x + eg) x - br 0 g , (14) 

r 2 g = r 0 x + (r 0 x + eg) g , (15) 

where 

r 2 = x 2 + bg 2 

Then f(t) and r( t ) are obtained from 

f = (x~ f Q g)/r 0 , 


r 


fr. 


g r ( 


For purposes of numerical control it may be advisable to let 
t - t 0 = hr , h=t 1 -t 0 , 


where [t 0 , t x ] is a time interval of interest. The corresponding interval for 
r is [O, 1] . Equations 12 and 13 still hold provided a 0 = r 0 , b 0 = b 2 = 0, 


3 


(t> -/i/r 0 ) h. 


a x - r 0 h, bj = h, a 2 - -/xh 2 /2r 0 2 , c - 

Equations 9 become 


x 



g " 


00 



i= 0 


which offer a further advantage when t - t x , i.e., when r = 1. 

Formulas for the radii of convergence of the f and g series are given by 
Moulton. 
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